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Interest  in  thermoelectrics  for  waste-heat  recovery  and  localized  cooling  has  flourished  in  recent  years,  but 
questions  about  cost  and  scalability  remain  unanswered.  This  work  investigates  the  fabrication  costs  and 
coupled  thermal  and  electrical  transport  factors  that  govern  device  efficiency  and  commercial  feasibility  of  the 
most  promising  thermoelectric  materials.  For  30  bulk  and  thin  film  thermoelectric  materials,  we  quantify  the 
tradeoff  between  efficiency  and  cost  considering  electrical  and  thermal  transport  at  the  system  level,  raw 
material  prices,  system  component  costs,  and  estimated  manufacturing  costs.  This  work  neglects  the  cost  of 
heat,  as  appropriate  for  most  waste-heat  recovery  applications,  and  applies  a  power  generation  cost  metric  in 
$/W  and  a  cooling  operating  cost  metric  in  $/l<Wh.  The  results  indicate  material  costs  are  too  high  for  typical 
thermoelectric  power  generation  applications  at  mean  temperatures  below  135  °C.  Above  275  C,  many  bulk 
thermoelectric  materials  can  achieve  costs  below  $1/W.  The  major  barrier  to  economical  thermoelectric  power 
generation  at  these  higher  temperatures  results  from  system  costs  for  heat  exchangers  and  ceramic  plates.  For 
cooling  applications,  we  find  that  several  thermoelectric  materials  can  be  cost  competitive  and  commercially 
promising. 
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Nomenclature 

LCOE 

levelized  cost  of  electricity  ($/kWh) 

71 

total  plate  area  (m2) 

m 

max 

ratio  of  load  resistance  to  device  resistance 
subscript  designating  maximum  value 

Act 

thermoelectric  leg  i  cross  sectional  area  (m2) 

n 

subscript  designating  n-type  semiconductor 

Ate 

total  thermoelectric  material  cross  sectional  area  (m2) 

opt 

subscript  designating  optimum  value 

C 

total  capital  cost  ($) 

p 

subscript  designating  p-type  semiconductor 

C" 

cost  by  volume  ($/m3) 

Pgen 

generator  electric  power  output  (W) 

C" 

cost  by  area  ($/ m2) 

Psup 

electrical  power  supplied  (W) 

CB 

bulk  raw  material  cost  ($/ kg) 

Qc 

heat  transferred  through  cold  side  heat  exchanger  (W) 

Ce 

price  of  electricity  ($/kWh) 

Qh 

heat  transferred  through  hot  side  heat  exchanger  (W) 

Chx 

heat  exchanger  cost  ($/(W/K)) 

Q.TE 

heat  transferred  through  legs  (W) 

Cm.b 

manufacturing  cost  to  process  bulk  material  ($/kg) 

Oil 

heat  transferred  between  legs  (W) 

Cma 

manufacturing  cost  to  process  material  by  area  ($/ m2) 

r 

amortization  (discount)  rate  (h_1) 

COP 

coefficient  of  performance 

R 

electrical  resistance  (Q) 

F 

fill  factor 

Rl 

electrical  load  resistance  (£2) 

G 

generator  cost-performance  metric  (S/W) 

S 

Seebeck  coefficient  (V/K) 

Grain 

minimum  cost  for  thermoelectric  device  ($/W) 

Sp„ 

difference  in  Seebeck  coefficients  of  p-  and  n-type  legs 

Go 

cost-performance  scaling  metric  for  thermoelectric 

(V/K) 

device  ($/W) 

T 

temperature  (K) 

H 

cooler  operating  cost  ($/kWh) 

Tc 

cold  side  temperature  (K) 

I 

electrical  current  (A) 

Th 

hot  side  temperature  (K) 

k 

thermal  conductivity  (W/m-K) 

Tm 

mean  junction  temperature  (K) 

I<c 

thermal  conductance  of  cold  side  heat  exchanger  (W/ 

Ti 

hot  junction  temperature  (K) 

K) 

t2 

cold  junction  temperature  (K) 

I<h 

thermal  conductance  of  hot  side  heat  exchanger  (W/I<) 

AT 

junction  temperature  difference  (K) 

K, 

total  thermal  conductance  (W/K) 

U 

heat  transfer  coefficient  of  heat  exchanger  (W/m2  K) 

Kje 

thermal  conductance  of  thermoelectric  legs  (W/I<) 

XT 

module  or  material  thermoelectric  figure-of-merit 

I<» 

thermal  conductance  of  parallel  path  (W/I<) 

>1 

efficiency 

L 

leg  length  (m) 

P 

density  (kg/m3) 

Ldr 

diminishing  returns  optimum  leg  length  (m) 

G 

electrical  conductivity  (S/m) 

lt 

ratio  of  thermal  conductivity  to  heat  transfer  coeffi¬ 
cient  of  heat  exchanger  (m) 

1.  Introduction 

Thermoelectric  devices  are  used  in  power  generation  and 
cooling  applications  to  either  convert  heat  into  electricity  or  to 
pump  heat.  To  date,  thermoelectric  technology  has  been  con¬ 
strained  to  applications  that  include  vehicle  waste-heat  recovery 
prototypes,  space  vehicle  power  sources,  seat  coolers,  solid-state 
refrigerators,  and  temperature  control  in  laboratory  equipment. 
Thermoelectric  conversion  has  received  renewed  interest  due  to 
the  development  of  better-performing  materials  and  their  poten¬ 
tial  to  improve  the  efficiency  of  combustion  systems  through 
waste-heat  recovery  [1-3].  Improvements  in  material  performance 
are  ongoing  [4-6]  since  conversion  efficiencies  of  typical  thermo¬ 
electric  materials  remain  below  10%  [3,7], 

Thermoelectric  devices  are  attractive  because  they  have  no 
moving  parts  and  can  be  silent,  reliable,  and  versatile.  Since 
multiple  thermoelectric  n-  and  p-type  couples  can  be  connected 
in  series,  a  thermoelectric  device  can  be  smaller  than  a  computer 
chip  or  larger  than  a  solar  panel.  Nevertheless,  considerable 
technical  challenges  remain.  Existing  thermoelectric  converters 
have  lower  conversion  efficiencies  than  alternatives  such  as 
organic  Rankine  cycles.  While  some  applications  are  well  estab¬ 
lished,  in  particular  power  sources  for  space  vehicles,  a  variety  of 
terrestrial  applications  -  such  as  high  temperature  waste-heat 
recovery  -  are  yielding  new  concerns  about  reliability  and  dur¬ 
ability.  Challenges  related  to  sublimation,  oxidation,  thermal  and 
electrical  interface  degradation,  and  thermal  expansion  mismatch 
remain  critical  for  applications  of  thermoelectric  devices  [8-15], 

Thermoelectric  technologies  face  several  additional  commer¬ 
cialization  challenges.  The  cost  of  many  thermoelectric  materials 


may  be  prohibitively  high,  largely  due  to  the  use  of  tellurium  and 
germanium  in  the  most  common  contemporary  thermoelectric 
applications  [16],  Safety  concerns  preclude  toxic  materials  such  as 
lead  [17],  The  weight  and  specific  power  of  thermoelectric  devices 
are  particularly  important  for  mobile  applications  such  as  vehicle 
waste-heat  recovery  [18],  yet  few  strategies  for  weight  reduction 
have  been  proposed.  While  government  funding  for  thermoelec¬ 
tric  technology  research  and  development  has  expanded  signifi¬ 
cantly  in  the  past  decade,  the  relative  lack  of  private  sector 
familiarity  with  the  technology  made  early  stage  financing  for 
companies  comparatively  slow  to  follow  [19].  There  has  been  a 
rapid  development  of  materials  promising  higher  efficiencies,  in 
many  cases  through  the  use  of  nanostructuring  and  novel  crystal 
structures  [1,4,6],  Recent  prototypes  demonstrated  the  feasibility 
of  thermoelectric  devices  for  widespread,  terrestrial,  scalable 
applications  [20-23],  Multiple  start-up  companies  have  received 
recognition  and  funding  while  also  forming  partnerships  with 
established  academic  and  industrial  research  teams  [15,18,24-26], 
The  continued  development  and  deployment  of  scalable  thermo¬ 
electric  devices  depend  on  the  device  cost  and  energy  conversion 
performance  [1,27], 

There  is  an  urgent  need  for  a  comprehensive  assessment  of 
the  commercial  feasibility  of  the  most  promising  thermoelectric 
materials  while  also  considering  relevant  manufacturing  and 
system  costs.  Cost-performance  analyses  of  this  type  have  been 
conducted  for  other  energy  technologies  such  as  batteries  and 
photovoltaics  by  considering  material  extraction  costs  [28,29],  For 
thermoelectric  devices,  there  have  been  multiple  evaluations  of 
the  coupled  thermal  and  electrical  module  performance  without 
considering  cost  [30,31],  The  importance  of  cost  was  considered  by 
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normalizing  the  material  figure  of  merit  by  the  raw  material  cost, 
and  the  costs  for  specialized  applications  such  as  a  marine  waste 
incinerator  have  also  been  determined  [17,27,32,33],  A  method  for 
determining  cost  using  the  thermal  energy  cost  and  the  module 
construction  cost  for  specified  material  thicknesses  has  been 
presented  [32],  A  recent  investigation  of  the  cost-efficiency  trade¬ 
offs  of  thermoelectric  generators  first  determined  the  geometry 
that  optimizes  the  thermal  and  electrical  performance  of  the 


device  and  then  calculated  the  material  cost  associated  with  the 
optimized  dimensions  [27],  The  work  highlighted  the  importance 
of  optimizing  modules  in  conjunction  with  the  heat  source  and 
sinks  in  the  overall  system.  An  important  extension  to  the 
optimization  approach  is  made  by  evaluating  the  coupled  thermal, 
electrical,  and  economic  factors  simultaneously  [34],  This  method 
captures  the  ability  of  a  moderately  efficient  device  with  low 
capital  cost  to  perform  better  in  the  market  than  a  highly  efficient 


Table  1 

Material  identification  table.  Material  type  is  indicated  in  the  first  column  and  by  material  name  color  coding:  red  for  chalcogenides  and  SiGe,  blue  for  silicides,  purple  for 
clathrates,  green  for  skutterudites,  black  for  oxides,  yellow  for  half-Heuslers,  and  grey  for  others.  (For  interpretation  of  the  references  to  color  in  this  table  caption,  the  reader 
is  referred  to  the  web  version  of  this  paper.) 


Material  type 

ID 

Material  name 

Manufacturing 

type 

Material  cost 

(S/kg) 

Evaluated  for  low  or  high  temp. 

scenario 

Ref. 

ZTm 

material 

Fopt 

f-opt 

(mm) 

Scenario 

Temp. 

Chalcogenides  and 

1 

Bi2Te3 

Bulk 

110 

0.74 

0.18 

4.53 

Low 

135,36] 

SiGe 

2 

Bio.52Sbi.4sTe3 

Bulk 

125 

1.05 

0.21 

4.41 

Low 

]37] 

3 

Bio.52Sbi.4sTe3 

Nanobulk 

125 

1.52 

0.29 

3.47 

Low 

[37] 

4 

Bio.54Teo.46 

Nanowire 

84 

0.02 

0.07 

1.10 

Low 

[38] 

5 

( Nao.0283Pbo.g45Teo.9733 ) 
(Ag1.nTeo.555) 

Nanobulk 

81 

1.45 

0.34 

3.01 

High 

[39] 

6 

Bi-doped  PbSeo.9sTeo.02/PbTe 

Superlattice 

55 

1.96 

0.02 

0.31 

Low 

[40] 

7 

AgPb,8SbTe2o 

Bulk 

84 

1.31 

0.26 

3.59 

High 

[41] 

8 

SiGe 

Bulk 

679 

0.30 

0.07 

2.66 

High 

[42] 

9 

SisoGe2o 

Nanobulk 

371 

0.53 

0.13 

3.39 

High 

[42] 

10 

SiGe 

Nanowire 

679 

0.22 

0.06 

1.59 

Low 

[43] 

Silicides 

11 

Mg2Sio.85Bio.15 

Nanobulk 

6.67 

0.67 

0.74 

29.5 

High 

[44] 

12 

Mg2Sio.6Sno.4 

Bulk 

4.04 

1.05 

~1 

45.5 

High 

[45] 

13 

Si 

Nanobulk 

3.09 

0.21 

0.71 

94.5 

High 

[46] 

14 

Si 

Nanowire 

3.09 

0.72 

0.09 

3.38 

Low 

[47] 

15 

MnSii.75 

Bulk 

1.46 

0.05 

~1 

37.5 

Low 

[48] 

16 

Mni5Si28 

Nanobulk 

1.51 

0.07 

~1 

35.7 

Low 

[48] 

Clathrates 

17 

Ba8Gai6Ge28Zn2 

Bulk 

615 

0.48 

0.13 

1.50 

High 

[49] 

18 

Ba8Gai6Ge3o 

Bulk 

644 

0.36 

0.11 

1.65 

High 

[50] 

19 

Ba7SriAli6Si3o 

Bulk 

1.64 

0.09 

~1 

38.6 

High 

[51] 

Skutterudites 

20 

CeFe4Sbi2 

Bulk 

37 

0.77 

0.28 

8.34 

High 

[52] 

21 

Yb0.2lno.2Co4Sbi2 

Bulk 

24 

0.93 

0.31 

10.6 

High 

[53] 

22 

Cao.i8C03.97Nio.03Sb12.40 

Bulk 

13 

0.77 

0.39 

17.6 

High 

[54] 

Oxides 

23 

(Zn0.98Alo.o2)0 

Bulk 

2.30 

0.08 

0.48 

58.6 

High 

[55] 

24 

Ca2.4Bio.3Na0.3Co409 

Bulk 

30 

0.13 

0.43 

7.17 

High 

[56] 

25 

InGaZnO 

Nanowire 

511 

0.07 

0.04 

1.59 

Low 

[39] 

26 

Nao.7Co02-8 

Bulk 

36 

0.52 

0.22 

12.7 

High 

[57] 

Half  Heuslers 

27 

Zr0.25Plf0.25Ti0.5NiSn0.994Sb0.006 

Bulk 

9.71 

1.38 

0.49 

19.6 

High 

[58] 

28 

Zro.5Hfo.5Nio.8Pdo.2Sno.99Sbo.01 

Bulk 

8.51 

0.69 

0.46 

15.4 

High 

[59] 

29 

Tio.sHf0.2NiSn 

Bulk 

10.70 

0.41 

0.41 

16.5 

High 

[60] 

Other 

30 

PEDOT:PSS 

Polymer 

0.34 

0.01 

~8 

101 

Low 

[61] 

*See  Table  2  below  for  scenario  values. 

1 :  Materials  with  Fopt~  have  high  thermal  conductivities  which  suggests  that  heat  exchangers  with  larger  (/-values  should  be  used. 

8:  The  cost  of  polymers  and  their  manufacturing  technique  are  so  different  from  other  thermoelectric  materials  that  a  different  geometry  than  conventional  TE  should  be 
considered. 


Table  2 

Scenario  parameters,  device  temperatures,  (/-values,  and  heat  exchanger  (HX)  costs.  Representative  applications  are  provided  for  each  temperature  scenario. 


Applications 

Scenario  temperature 

Thermoelectric 

cooling 

(Cooling) 

Low  temperature 
recovery 

Scenario  #1 

Low 

Solar  thermal 
generator 

Scenario  #2 

Medium-Low 

Automotive  exhaust  heat 
recovery 

Scenario  #3 

Medium-High 

Industrial  furnace  heat 
recovery 

Scenario  #4 

High 

Hot  side  temperature, 

Th  (°C) 

15 

100 

250 

500 

800 

Cold  side 
temperature, 

Tc  CC) 

5 

20 

20 

50 

50 

Average  temperature, 

Tm  (°C) 

10 

60 

135 

275 

425 

Hot  side  (/-value, 

(/„  (W/m2  K) 

00 

102 

102 

120 

120 

Cold  side  (/-value, 

Uc  (W/m2  K) 

00 

105 

105 

105 

105 

HX  &  Plate  Costs 

($/(  W/K)) 

18.48 

18.48 

18.48 

18.48 
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device  with  higher  capital  cost.  While  this  prior  work  made 
progress  on  coupling  efficiency  and  cost,  the  actual  values  and 
thus  large  impacts  of  manufacturing  costs  and  the  costs  of  key 
system  components  have  not  been  considered.  This  is  particularly 
important  for  the  newer  bulk  and  thin  film  materials  for  which 
manufacturing  costs  can  be  high,  and  new  system  designs  might 
be  required  to  optimize  heat  exchange. 

The  present  work  provides  a  detailed  investigation  of  the 
materials,  manufacturing,  and  system  costs  for  30  leading  thermo¬ 
electric  materials  [35-61]  considering  the  conversion  efficiency. 
We  combine  these  using  cost  metrics  that  account  for  the  coupled 
cost  and  system  design  implications  and  trade-offs  [34].  These 
metrics  enable  an  assessment  of  which  materials  are  most  pro¬ 
mising  for  heat  recovery  and  cooling  and  incorporate  an  optimiza¬ 
tion  of  material  thickness  and  fill  factor.  The  present  work  sets 
targets  for  thermoelectric  cooling  and  power  generation.  The  cost 
of  power  generation,  in  $/W,  is  reported  by  combining  material 
properties,  device  physics,  and  raw  material  costs,  all  of  which  are 
reported  for  seven  material  classes  presented  in  Table  1  and  in 
Appendix  F.  The  analysis  is  performed  for  five  operating  tempera¬ 
tures  shown  in  Table  2  to  reflect  the  myriad  of  potential  thermo¬ 
electric  applications.  The  difference  in  manufacturing  costs  for 
bulk  and  thin  film  materials  influences  the  decision  on  whether  to 
use  novel  materials.  A  comparison  to  existing,  competitive  tech¬ 
nologies  highlights  the  cost  of  current  state-of-the-art  thermo¬ 
electric  technologies  and  device  performance. 


2.  Methods 

This  work  determines  the  appropriate  device  geometry  based 
on  cost  and  performance.  The  fill  factor,  F,  and  the  leg  length,  L,  of 


a  thermoelectric  device  are  the  dominant  design  parameters.  The 
fill  factor  is  the  ratio  of  the  area  covered  by  the  active  thermo¬ 
electric  material  to  the  plate  area,  A,  as  illustrated  in  Fig.  1.  While  it 
goes  against  traditional  thermoelectric  module  architecture,  it  is 
possible  to  envision  devices  with  F>  1  (see  Section  3.1).  The 
inverse  of  the  fill  factor  is  known  as  the  thermal  concentration 
factor  [62].  By  reducing  the  fill  factor,  less  active  material  is  used, 
and  the  overall  cost  of  the  device  can  be  reduced.  In  many 
applications,  the  fill  factor  is  fixed  to  be  above  some  minimum 
by  other  system  constraints  (e.g.,  mechanical  robustness,  voltage/ 
areal  requirements).  By  reducing  the  leg  length  at  this  fixed  fill 
factor,  less  material  is  used,  and  the  cost  decreases.  However,  for  a 
fixed  heat  exchanger  and  fixed  F,  as  the  thermoelectric  leg  length 
is  reduced  below  some  optimum,  the  temperature  drop  across  the 
device  decreases,  and  the  output  power  decreases.  This  competi¬ 
tion  between  power  output  and  material  volume  suggests  the 
existence  of  an  optimum  geometry  which  minimizes  cost  and 
maximizes  power  output.  Therefore,  the  most  appropriate  metric 
for  thermoelectric  generators  should  be  expressed  on  a  $/W  basis 
as 


where  Pgen  is  the  electric  power  output  of  the  generator,  and  C  is 
the  capital  cost  of  the  device.  The  capital  cost  includes  the 
material,  manufacturing  and  heat  exchanger  costs  which  are 
described  in  detail  in  Section  2.2.  Cost  of  the  heat  available  for 
conversion  is  excluded.  The  contribution  of  heat  exchanger  and 
insulating  ceramic  plate  costs  are  considered  in  Section  3.2.  Both  C 
and  Pgen  are  strongly  dependent  on  L  and  F.  A  cost-performance 
optimization  should  be  performed  on  C  rather  than  on  C  or  Pgen 
individually  [34]. 


Thermoelectric  area,  Aj-e  =  Acn  +  Ac 
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Fig.  1.  Diagram  of  thermoelectric  module  and  the  equivalent  thermal  circuit.  The  n-  and  p-type  thermoelectric  materials  are  connected  electrically  in  series  with  a  metal 
(shown  in  orange).  The  thermoelectric  legs  are  thermally  in  parallel  between  electrically  insulating  ceramic  plates  (shown  in  light  grey).  The  analysis  is  conducted  for 
constant  hot-  and  cold-side  temperatures  TH  and  Tc.  There  are  two  thermal  conductance  terms  between  the  junction  temperatures  T i  and  T2  representing  heat  transfer 
through  the  thermoelectric  material  legs  and  in  the  gap  between  the  legs.  The  internal  heat  generation  term,  I2R/2,  is  the  Joule  heating  component;  the  Peltier  term  is  ISpnT. 
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A  similar  argument  can  be  made  for  thermoelectric  coolers. 
In  order  to  reduce  the  capital  cost  of  the  module,  shorter  leg 
lengths  are  used.  However,  this  reduces  the  efficiency  of  the 
thermoelectric  cooler.  Since  thermoelectric  coolers  consume  elec¬ 
tricity  and  provide  a  cooling  service,  the  cost  of  electricity  is 
important,  so  a  thermoelectric  module  with  a  large  coefficient  of 
performance,  COP,  is  desirable.  The  COP  is  the  ratio  of  the  cooling 
power  removed  from  the  cold  side  to  the  electric  power  con¬ 
sumed.  Electricity  is  purchased  on  a  $/l<W  he  basis,  so  the  most 
appropriate  metric  is  the  operating  cost  expressed  in  units  of 
S/kW  hth  where  the  subscripts  represent  electric  and  thermal 
(cooling)  energy,  respectively.  The  capital  cost  of  the  thermo¬ 
electric  cooler  is  amortized  over  the  lifetime  and  added  to  the 
operating  cost.  Operating  continuously  and  purchasing  electricity 
hourly  over  the  lifetime  allows  for  a  simplification  in  the  amorti¬ 
zation  where  the  capital  cost  is  multiplied  by  the  amortization 
(discount)  rate,  r,  on  an  hourly  basis  (see  Appendix  B).  This  allows 
for  both  capital  cost  and  operating  costs  to  be  considered  for 
cooling  applications.  The  appropriate  operating  cost  metric  is 

H  =  COP+rCOPPsup’  (2) 


temperatures  and  working  fluids).  The  [/-value  of  the  heat 
exchangers  and  the  thermal  conductance  are  related  to  each  other 
by  the  heat  transfer  area,  I<=UA,  and  the  values  [68]  used  in  this 
analysis  along  with  application  reservoir  temperatures  are  sum¬ 
marized  in  Table  2.  For  the  scenarios  selected  in  this  analysis,  the 
appropriate  hot  side  and  cold  side  exchangers  have  approximately 
the  same  [/-value,  so  matched  heat  exchanger  simplifications  (i.e., 
Kh=Kc=I<)  are  appropriate.  However,  for  exactness  these  simpli¬ 
fications  are  not  used  in  this  analysis.  Ultimately,  these  series 
conductances  relate  the  hot-  and  cold-side  reservoir  temperatures, 
Th  and  Tc,  to  the  hot-side  and  cold-side  junction  temperatures, 
7j  and  T2. 


2.1.1.  Power  generation 

The  temperature  difference  governs  the  energy  that  can  be 
converted  into  a  thermoelectric  generator.  For  the  electrical  load 
matched  condition,  energy  balances  of  the  device  at  the  hot  and 
cold  junctions  yield  a  system  of  coupled  non-linear  equations: 


Snn(7J  — T2)T\ 

KH(TH-Tr)-KT(Tr-T2)--pn 


2  R 


S2pn(Tt-T2)2 

M 


=  0 


(3a) 


where  Ce  represents  the  price  of  electricity,  and  Psup  is  the  electric 
power  supplied  to  the  device  to  maintain  a  desired  temperature 
difference  across  the  module.  The  price  of  electricity  varies.  In  the 
United  States,  the  average  industrial  user  purchases  electricity  at  a 
price  of  $0.068/kW  hc,  and  residential  users  purchase  electricity  at 
a  price  of  $0.12/1<W  he.  The  average  price  of  electricity  is  taken  as 
$0.098/kW  hc  in  this  analysis  and  is  bounded  by  the  industrial  and 
residential  prices  (i.e.,  the  error  bars).  These  values  are  the  average 
U.S.  electricity  prices  reported  by  the  Energy  Information  Agency 
at  the  time  of  article  submission.  The  product  of  the  COP  and  Psup  is 
the  thermal  cooling  power,  in  l<Wth,  and  yields  the  thermal  energy 
removed  in  kW  hth  once  divided  by  r. 

2.1.  Device  physics 

Both  of  these  cost  metrics  rely  on  first  modeling  the  physics  of 
the  thermoelectric  device.  Fig.  1  depicts  a  dual-leg  thermoelectric 
module  with  the  legs  connected  electrically  in  series  and  ther¬ 
mally  in  parallel.  The  electrical  and  thermal  transport  for  this 
geometry  has  routinely  been  solved  [30,63-66],  The  total  thermal 
conductance  of  the  thermoelectric  module,  I(T,  is  the  sum  of  the 
thermal  conductance  through  the  active  material,  KTE,  and  a 
parallel  thermal  leakage  conductance,  K„.  The  heat  transferred 
between  the  legs,  Q,|,  is  described  by  which  is  the  sum  of 
conduction,  convection,  and  radiation  contributions.  For  most 
applications,  conduction  is  the  dominant  heat  transfer  mechanism, 
but  radiation  and  convection  can  be  significant.  As  a  result  of  these 
three  contributions,  K„  and  I<T  are  non-linear  functions  of  leg 
length  L.  A  description  of  how  kj,  is  handled  in  this  analysis  can 
be  found  in  Appendix  D.  The  heat  flowing  through  the  thermo¬ 
electric  legs,  QT£,  is  augmented  by  the  Joule  heating  distributed 
throughout  the  thermoelectric  legs  and  Peltier  cooling/heating 
that  occur  locally  at  the  hot  and  cold  junctions. 

In  addition  to  the  thermal  conductance  of  the  module,  the 
thermal  conductance  of  the  hot-  and  cold-side  heat  exchangers,  I(H 
and  !<c,  respectively,  strongly  influence  the  overall  device  perfor¬ 
mance.  These  conductances  are  further  reduced  by  any  series 
temperature  drops  due  to  the  electrically  insulating  ceramic  plates 
and  module  thermal  contact  resistances.  The  effect  on  device 
performance  from  the  heat  exchangers  has  been  considered  [67], 
In  this  analysis,  the  contribution  from  the  metallization,  ceramic 
plates,  and  heat  exchangers  is  also  considered,  and  details  can  be 
found  in  Appendix  D.  The  values  of  the  thermal  conductances  are 
dependent  on  the  application  and  heat  exchangers  (i.e.,  reservoir 


Spn2(Tt-T2)2 
4  R 


-  Rh(Th~  T\)  +  Kc(T2  -Tc)  =  0, 


(3b) 


which  can  be  solved  numerically  for  the  exact  junction  tempera¬ 
tures.  In  this  expression,  Spn=Sp-Sn  is  the  difference  in  thermo¬ 
powers  of  the  two  legs,  and  R  is  the  internal  electrical  resistance  of 
the  module.  Under  the  approximations  that  the  electrical  contact 
resistance,  metallization,  and  shunt  resistances  are  negligible,  and 
the  thermoelectric  legs  have  the  same  electrical  conductivity 
<7p=<7n=er,  the  internal  electrical  resistance  is  R=4LI(crAF).  The 
analysis  is  conducted  assuming  the  properties  are  constant  within 
each  of  the  five  temperature  ranges  considered,  and  the  properties 
are  evaluated  at  the  mean  of  the  reservoir  temperatures  (i.e., 
7m=(r„+rc)/2).  Thermal  and  electrical  contact  resistances  are  not 
considered  as  these  vary  widely  and  are  unique  to  each  device; 
they  are  most  important  for  short  leg  lengths,  typically  <  1  mm. 
These  effects  can  significantly  reduce  the  performance  of  a  device, 
but  neglecting  these  terms  provides  a  lower  (optimistic)  bound  on 
the  $/W  cost  determined  for  each  material. 

The  amount  of  power  converted  by  a  thermoelectric  generator 
is 


Pgen  — 


SPn2(h-T2)2  m 

R  (m+1)2’ 


with 


(4) 


(5) 


as  the  ratio  of  the  load  resistance,  RL,  to  the  internal  electrical 
resistance,  R,  of  the  module.  The  thermoelectric  generator  can  be 
operated  at  maximum  power  when  m  —  1  (i.e.,  load  matched)  or  at 
maximum  efficiency  when  m  =  Vl  +ZTm.  It  has  been  recently 
reported  that,  due  to  thermal  feedback  from  the  load  resistance 
( e.g .,  increasing  the  load  resistance  increases  the  junction  tem¬ 
peratures  difference  slightly),  m  should  be  slightly  greater  than 
these  nominal  values  [31,69],  While  this  is  strictly  true,  the  effect 
is  a  higher  order  correction,  and  traditional  analysis  [30,31] 
suggests  that  one  of  these  m  conditions  should  be  selected.  The 
value  of  m  that  minimizes  the  $/W  is  m=  1  because  this  condition 
maximizes  the  denominator  of  Eq.  (1).  Therefore,  the  load  resis¬ 
tance  is  set  to  exactly  equal  the  internal  resistance  [34], 

In  this  work,  ZT  is  the  module  figure-of-merit 


7T  =  ^ElL 
~  KtR 


(6) 
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and  Tm=(T^  +  T2)l2  is  the  average  of  the  junction  temperatures. 
Eq.  (6)  is  similar  to  the  material  figure  of  merit  ZT=S2crTlk,  with 
the  material's  thermopower,  S,  electrical  conductivity,  a,  and 
thermal  conductivity,  k.  The  module  and  material  ZT  are  related 
as  discussed  in  Appendix  A.  For  this  manuscript,  the  material  ZT 
will  be  used  to  identify  the  materials  in  the  figures  while  the 
module  ZT  is  used  in  the  analysis  under  the  simplification  that  the 
n-  and  p-  type  legs  have  the  same  thermopower  magnitude,  same 
electrical  conductivity,  and  same  thermal  conductivity. 

Using  the  electric  power  and  the  cost  of  the  device,  Eq.  (1)  can 
be  evaluated  in  units  of  $/W  as  a  function  of  L  for  a  chosen  F  as 
shown  in  Fig.  2a  for  a  doped  bismuth  telluride  example  (Material 
ID  #2).  By  selecting  a  fill  factor,  an  optimum  length,  Lopt,  exists 
which  minimizes  G.  The  heat  flowing  into  the  hot  side,  Qh,  and  the 
thermal  efficiency,  >]=PgenlQH,  are  determined  as  a  function  of  leg 
length  as  shown  in  Fig.  2b.  Using  Eqs.  (3)— (5),  the  electric  power 
generated  is  determined  as  a  function  of  leg  length  as  shown  in 


Length  (m)  Length  (m) 


Fig.  2.  Example  of  optimum  leg  length  Lop,  and  cost  C  minimization  for  power 
generation  using  Material  ID  #2  in  Table  1  (bulk  Bio.52Sbi  4sTe3)  with  Fopt= 0.21  and 
parameters  in  the  low  temperature  scenario  of  Table  2.  The  leg  length  that  results 
in  the  minimum  cost  (a,  Lopt=4.4  mm)  is  not  the  same  length  that  maximizes  the 
efficiency  (b,  Lmax,p=24  mm)  or  output  power  (c,  LmaxP—7.3  mm).  Note  (a)  is  on  a 
linear  scale  while  (b)  and  (c)  are  on  a  log  scale.  The  black  and  blue  n  curves  in 
(b)  are  nearly  indistinguishable.  (For  interpretation  of  the  references  to  color  in  this 
figure  caption,  the  reader  is  referred  to  the  web  version  of  this  paper.) 


Fig.  2c.  The  leg  length  which  minimizes  G  is  different  from  the  leg 
length  that  maximizes  efficiency  or  the  leg  length  that  maximizes 
the  electric  power  output.  Fig.  2  shows  the  difference  between 
designs  optimized  for  efficiency  or  power  and  an  alternative 
design  that  optimizes  the  ratio  of  cost  to  power  output. 

The  cost,  efficiency,  and  electric  power  output  are  functions  of 
the  load  condition  m  as  discussed  in  Appendix  A  and  as  illustrated 
by  the  different  color  lines  in  Fig.  2.  In  situations  where  there  is 
excess  heat  (i.e.,  fuel/heat  is  free)  the  maximum  power  condition 
m  =  1  is  most  appropriate  and  yields  the  lowest  cost.  In  situations 
where  heat  is  costly,  the  maximum  efficiency  condition  may  be 
more  appropriate.  Evaluating  fuel  costs  is  beyond  the  scope  of  this 
work  but  has  been  investigated  elsewhere  [32],  The  unmatched 
condition  (i.e.,  m  is  arbitrary)  may  also  be  acceptable  depending  on 
the  application  or  cost  of  power  electronic  conditioning  circuits 
necessary  to  load  match. 

Power  electronic  conditioning  circuits  match  the  resistance  of 
the  load  to  the  resistance  of  the  module  in  order  to  maintain  either 
the  maximum  power  condition  (m  =  1 )  or  the  maximum  efficiency 
condition  (m  =  +ZTm ).  The  cost-necessity  of  these  circuits  is 
evaluated  in  Fig.  3.  The  lack  of  a  proper  load  condition  results  in  an 
opportunity  cost  for  power  generation  which  is  plotted  as  a 
function  of  ZTm.  The  opportunity  cost  is  expressed  as  a  percentage 
of  the  capital  cost  of  the  device.  The  m  =  1  (horizontal  blue  line) 
case  is  the  ideal  situation  where  there  is  no  opportunity  cost. 
Other  arbitrary  load  conditions  (horizontal  red  lines)  represent 
different  degrees  of  opportunity  costs.  The  black  curve  in  Fig.  3  is 
the  percent  difference  between  G  with  m  =  Vl  +ZTm  and  G  with 
ni=l;  this  is  the  opportunity  cost  for  operating  at  the  maximum 
efficiency  condition  rather  than  the  load  matched  condition.  The 
black  circles  represent  the  position  of  materials  considered  in  this 
study  along  the  maximum  efficiency  line.  For  most  materials  in 
this  study,  ZTm  is  low  enough  that  operating  at  the  maximum 
efficiency  or  maximum  power  makes  little  difference.  However, 
at  larger  ZTm  this  difference  becomes  significant.  For  a  device 
operating  at  a  condition  other  than  m  =  l,  the  opportunity  cost  is 
higher,  and  the  added  cost  of  power  electronics  should  be  consid¬ 
ered.  For  example,  if  a  material  has  a  ZTm=  1,  and  throughout  normal 


Fig.  3.  Cost  associated  with  non-optimized  load  conditions.  The  opportunity  costs 
are  percentages  of  the  total  capital  cost  of  the  module.  Designing  the  module  to 
operate  at  the  maximum  power  condition  or  introducing  a  power  electronics 
conditioning  circuit  to  load  match  has  appreciable  cost  benefit.  The  horizontal  lines 
represent  specific  load  conditions  m=RLIR:  (blue,  m  =  1)  maximum  power  condi¬ 
tion  (black,  m  =  VI  +ZTm),  maximum  efficiency  condition  (red,  m  =  1.50, 1.75,  2.00) 
unmatched  load  conditions.  The  circles  represent  the  practical  maximum  efficiency 
conditions  of  the  materials  investigated  in  this  study.  The  optimum  load  is  the 
matched  load  (m=l),  which  produced  the  maximum  power  and  the  zero 
opportunity  cost.  (For  interpretation  of  the  references  to  color  in  this  figure 
caption,  the  reader  is  referred  to  the  web  version  of  this  paper.) 
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operation  m  fluctuates  to  1.5  for  a  large  fraction  of  generating  time, 
then  an  opportunity  cost  of  4%  of  the  capital  cost  is  incurred  in  the 
absences  of  conditioning  circuits.  If  the  conditioning  circuit  costs  less 
than  4%  of  the  capital  cost  of  the  device,  then  the  additional  power 
delivered  warrants  the  cost  of  the  conditioning  circuit.  This  varies 
with  each  material/system  as  they  have  different  capital  costs.  For 
inexpensive  systems  the  additional  system  cost  of  conditioning 
circuits  is  not  justified  whereas  this  cost  may  be  justified  for 
expensive  systems. 

2.1.2.  Levelized  cost  of  energy  considerations 

The  levelized  cost  of  energy  (LCOE)  is  often  used  to  compare 
the  value  of  different  energy  generating  technologies;  a  lower 
LCOE  value  corresponds  to  a  lower  cost  to  generate  electricity  [70], 
This  metric  is  a  function  of  several  variables  including  the  over¬ 
night  capital  cost,  the  ongoing  fuel  costs,  and  the  recurring 
maintenance  costs.  These  costs  are  spread  over  the  lifetime  of 
the  system  at  a  specified  capacity  factor  (uptime)  and  financial 
discount  rate.  The  units  of  LCOE  are  typically  $/kW  h,  and 
determining  the  LCOE  value  is  application-specific  because  of  the 
number  (and  specificity)  of  input  parameters.  The  analysis  in  this 
manuscript  only  calculates  the  overnight  capital  cost  of  the 
generator  (related  to  G),  which  is  only  one  input  parameter  for 
the  LCOE  value.  A  simplified  LCOE  method  is  presented  in 
Appendix  C  and  provides  a  way  of  converting  the  $/W  value  into 
a  S/kW  h  value  using  more  specific  parameters  which  are  beyond 
the  scope  of  this  manuscript. 


2.1.3.  Thermoelectric  cooling 

The  device  physics  of  thermoelectric  coolers  are  similar  to 
those  of  generators  except  that  electric  power  is  supplied  to 
support  a  temperature  difference  rather  than  being  generated. 
It  is  usually  desirable  to  operate  thermoelectric  coolers  at  the 
maximum  COP  since  electricity  is  continually  purchased.  For 
thermoelectric  cooling  the  goal  is  to  maintain  a  specified  tem¬ 
perature  difference  at  the  junctions  of  the  module.  The  hot-  and 
cold-side  heat  exchangers  act  as  additional  thermal  resistances. 
The  cold-side  junction  T2  needs  to  be  colder  than  7C,  and  the  hot 
side  junction  7,  needs  to  be  hotter  than  TH.  This  over  cooling/ 
heating  places  an  extra  burden  on  the  thermoelectric  module 
design  as  the  cold  side  junction  needs  to  be  cooled  below  the 
target  reservoir  temperature.  Removing  this  burden  for  simplifica¬ 
tion,  ideal  hot  side  and  cold  side  heat  exchangers  are  considered 
here  (i.e.,  U=UH=Uc=oo).  However,  the  series  thermal  resistances 
of  the  ceramic  plates  and  metal  shunts  are  still  included,  so  1<H 
and  Kc  are  finite  (see  Appendix  D).  This  simplification  fixes 
the  temperature  of  the  exteriors  of  the  ceramic  plates  to  be  the 
reservoir  temperatures.  The  fill  factor  F  is  approximately  one 
because  the  device  should  be  covered  in  as  much  active  thermo¬ 
electric  material  as  possible  in  order  to  pump  heat. 

Under  these  simplifications  ( U=oo ,  F=l),  the  mean  junction 
temperature  is  approximated  as  the  mean  reservoir  temperature, 
7m «  (Th+Tc)/2,  enabling  an  explicit  solution  for  the  junction 
temperatures  by  performing  an  energy  balance  on  the  cold  side 
junction.  This  yields  a  quadratic  equation: 


(  S2pnVT+ZTn  \  2 
V  KfVi+ZTro-iyy  2 

fSl„Tmy/ 1  +ZTm  + 1  Kc\ 

+[-Rurm^TF  K,+^) 
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for  which  the  smallest  root  is  the  cold  side  junction  temperature, 
and  the  hot  side  junction  temperature  is  determined  using 


Ti  ss  2Tm-T2.  These  junction  temperatures  are  used  to  determine 
the  COP  and  supplied  power  Psup: 


COP  = 


7,-72 


Vl+Z7m  — 7i/72 


(8) 


_  2Spn2(7,-72)2 
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Both  of  these  equations  are  functions  of  leg  length  and  can  be 
evaluated  numerically.  Finally,  using  these  equations,  the  operat¬ 
ing  cost  of  thermoelectric  cooling,  H  in  Eq.  (2),  can  be  evaluated. 


2.2.  Cost  considerations 


In  addition  to  the  device  physics,  both  material  and  manufac¬ 
turing  costs  influence  the  overall  cost  metric.  The  cost,  C  in  Eqs. 
(1 )  and  (2),  is  a  function  of  the  bulk  raw  material  cost  CB  ($/kg)  as 
summarized  in  Table  1.  The  manufacturing  cost  associated  with 
processing  bulk  material  CMB  ($/kg),  the  areal  manufacturing  cost 
Cma  ($/m2),  and  the  cost  of  both  heat  exchangers  Chx  ($/( W/K))  are 
found  in  Table  F.2  in  Appendix  F.  The  overall  cost  C  is  also 
dependent  on  the  heat  exchanger  heat  transfer  coefficient  U,  the 
density  of  the  active  material  p,  the  leg  length,  and  area  of  the 
ceramic  plates,  A: 


C  = 


(CB  +  CMB)pL+  cMA  \af+chxua 
V  c"  "  '~r~' ) 


(10) 


The  collection  of  terms  in  the  first  expression  is  a  volumetric  cost 
C"  ($/ m3),  and  the  second  term  is  an  areal  cost  C"  ($/m2),  both  of 
which  include  manufacturing  costs.  These  costs  are  determined 
based  on  a  conceptual  estimating  technique  [71]  since  these 
devices  are  largely  unprecedented.  This  technique  projects  costs 
for  future  products  which  do  not  exist  rather  than  gathering  prices 
for  components  of  a  device  already  in  production.  Rapid  iterations 
on  cost  estimates  completed  in  tandem  with  thermoelectric 
product  development  will  prove  valuable  as  devices  enter  the 
market. 

The  cost  of  the  extracted  raw,  bulk  material  CB  captures  the 
fundamental  differences  in  material  costs.  The  material  cost  is 
determined  from  the  2011  price  of  each  element  as  reported  by 
the  U.S.  Geological  Survey;  it  is  based  on  the  worldwide  produc¬ 
tion  of  the  element  and  represents  the  average  price  an  industrial 
consumer  would  pay  [72],  Thermoelectric  materials  are  typically 
composed  of  raw  materials  which  have  been  processed  to  a  purity 
level  of  99%  or  higher  [16].  Typical  costs  for  pure  materials  are 
presented  in  Appendix  F  and  can  be  significantly  higher  than  the 
raw  material  cost.  For  instance,  aluminum  is  an  abundant  material 
with  a  raw  material  cost  of  $2. 60/kg.  However,  the  cost  of  99.999% 
pure  aluminum  is  approximately  $300/kg  [73].  While  purification 
adds  to  the  material  cost  [16],  the  cost  of  purified  materials  is 
highly  uncertain.  It  can  vary  significantly  by  vendor  and  is  strongly 
influenced  by  rapidly  changing  market  and  processing  factors 
[14,59,74,75],  This  volatility  can  convolute  the  cost-performance 
analysis,  so  the  cost  to  purify  materials  is  omitted  here  without 
biasing  the  results  since  the  purification  conditions  apply  to  most 
minerals  [29],  A  comparison  between  pure  material  and  raw 
material  costs  was  conducted,  and  the  relative  cost  differences 
between  each  material  and  their  subsequent  ranking  did  not 
change  when  pure  material  costs  were  considered. 

Manufacturing  costs  are  divided  into  two  categories  based  on 
the  method  of  material  processing.  Processes  enacted  on  the  entire 
bulk  volume  of  material  such  as  ball  milling  and  spark  plasma 
sintering  have  costs  denoted  by  CMB.  Other  processes  like  dicing 
and  metallization  depend  on  the  area  of  material  processed;  these 
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costs  are  indicated  by  CMy \.  The  cost  is  calculated  by  dividing  the 
equipment  capital  cost  by  the  equipment  lifetime  and  throughput. 
The  exact  cost  and  throughput  values  are  available  in  Appendix  F. 
In  most  cases,  the  cost  and  throughput  of  equipment  were 
obtained  from  quotes  or  estimates  provided  by  equipment  ven¬ 
dors.  Throughput  is  highly  dependent  on  the  specific  application 
and  product,  so  it  is  not  easily  generalized  to  a  variety  of  materials 
and  applications.  An  activity  measure  [71]  has  been  assigned  to 
the  manufacturing  equipment  considered  here  which  assumes  the 
equipment  will  be  in  operation  year-round.  A  20-year  equipment 
lifetime  and  24  h  per  day,  365  days  per  year  operation  were 
assumed.  This  approach  considers  the  total  equipment  effective¬ 
ness  [76]  and  provides  the  lower  bound  on  manufacturing  cost. 
An  overall  equipment  effectiveness  which  incorporates  equipment 
shutdown  for  issues  like  repair  and  maintenance  could  also  be 
considered  [76].  The  effect  on  the  results  presented  here  are 
minimal  for  realistic  overall  equipment  effectiveness  values.  When 
the  analysis  is  conducted  assuming  equipment  operation  for  75% 
of  the  year,  the  resulting  increase  in  the  cost  metric  G  is  —5%  or 
lower  for  most  materials.  The  exception  is  the  superlattice  mate¬ 
rial  (Material  ID  #6)  for  which  the  increase  is  ~10%.  Additional 
costs  of  tooling,  maintenance,  repair,  assembly  (automated  or 
manual),  overhead,  and  labor  are  not  included  in  this  analysis. 

The  required  throughput  of  thermoelectric  material  in  a  man¬ 
ufacturing  process  will  depend  on  both  the  performance  of  the 
material  and  the  yield  of  the  process.  Fewer  legs  will  be  required 
for  higher  thermopower  materials,  and  less  material  must  be 
processed  as  yield  improves.  The  amount  of  active  material 
consumed  may  be  different  than  the  amount  of  active  material 
in  the  final  device  due  to  material  losses  during  manufacturing 
processes.  For  instance,  material  can  be  lost  during  dicing  pro¬ 
cesses,  deposited  on  chamber  walls  during  deposition,  or  removed 
through  etching.  The  impact  of  any  fractional  loss  of  the  material 
on  the  end  cost  would  be  the  same  as  any  changes  made  to  the 
equipment  effectiveness  as  described  above.  A  description  of  how 
to  adjust  the  manufacturing  cost  for  this  fractional  loss  is  provided 
in  the  Appendix  F.  In  order  to  provide  the  lower  bound  for  cost,  the 
results  reported  neglect  the  additional  cost  of  material  lost  due  to 
manufacturing  processes. 

Manufacturing  equipment  costs  are  approximate  and  estimated 
for  new  equipment.  Equipment  is  typically  tailored  specifically 
for  the  end-use  application,  so  communication  with  vendors  is 
required  for  accurate  quotes  and  precise  costs.  Particularly  in  early 
stage  research,  development,  and  production,  tools  built  in-house 
or  used/refurbished  equipment  may  be  used  to  reduce  costs.  In 
some  cases,  it  may  be  more  cost  effective  to  use  a  service  provider 
or  user  facility  for  certain  manufacturing  steps  rather  than 
purchasing  equipment  and  performing  the  process  in-house 
[36,77,78].  Equipment  and  manufacturing  costs  used  herein  sen/e 
as  points  of  reference;  real  manufacturing  costs  vary  by  material, 
process,  and  application. 

Table  1  summarizes  the  state-of-the-art  materials  with  pub¬ 
lished  and  peer-reviewed  thermoelectric  properties  used  in  the 
analysis.  They  are  from  the  following  classifications:  chalcogen- 
ides,  silicides,  clathrates,  half  Heusler  alloys,  skutterudites,  and 
oxides.  A  novel  polymer  material,  PEDOT:PSS,  has  recently  been 
proposed  for  thermoelectrics  and  is  also  included  here  [61],  The 
areal  manufacturing  costs  associated  with  the  polymer  material 
are  unique  since  a  polymer  device  would  likely  be  made  using  a 
process  like  screen  printing.  A  solution  printing  process  can  also  be 
used  for  the  metallization  layer  in  these  cases,  eliminating  the 
need  for  dicing.  Hence,  CMyt  for  polymer  or  solution-processed 
materials  could  be  significantly  lower  than  other  materials. 

Both  bulk  and  film  materials  are  considered  with  the  latter 
referring  to  thin  film,  superlattice,  and  nanowire  materials.  The 
raw  material  cost  for  each  thermoelectric  material  is  also  reported 


in  Table  1  along  with  the  references  for  material  properties  and 
characterization  results.  Materials  are  also  designated  by  their 
manufacturing  class:  bulk,  nanobulk,  nanowire,  superlattice,  and 
polymer.  Nanobulk  materials  are  ones  which  were  identified  by 
the  original  authors  as  having  nanoscale  grain  structures. 

The  limitations  of  manufacturing  techniques  in  turn  limit  the 
achievable  range  of  leg  lengths.  If  the  optimum  leg  length  based 
on  the  cost-performance  optimization  lies  outside  the  manufac¬ 
turing  limits,  entirely  new  manufacturing  methods  would  have  to 
be  developed  at  a  significant  additional  manufacturing  cost.  Bulk 
materials  should  be  limited  to  a  minimum  leg  length  of  500  pm 
as  defined  by  dicing  saw  limits  [79];  however,  after  executing 
the  analysis  none  of  the  bulk  materials  selected  for  this  study 
approached  this  dicing  limit.  Furthermore,  no  leg  length  limitation 
was  placed  on  nanowire  materials  since  nanowires  shorter  than 
500  pm  can  be  fabricated,  and  there  is  no  inherent  manufacturing 
limitation  to  making  long  thermoelectric  legs  with  nanowires.  The 
intrinsic,  single-nanowire  material  properties  reported  in  the 
references  considered  here  may  be  different  than  a  geometry  with 
multiple  nanowires,  but  those  considerations  are  case-specific  and 
beyond  the  scope  of  this  analysis. 

Potential  thermoelectric  power  generation  applications  span  a 
large  range  of  operating  temperatures.  Therefore,  the  analysis 
was  conducted  for  five  operating  temperatures  representing 
various  proposed  applications  as  shown  in  Table  2.  Temperature- 
dependent  material  properties  were  extracted  from  the  references 
listed  in  Table  1,  and  the  material  properties  at  the  mean  reservoir 
temperature  were  used  in  the  analysis. 


2.3.  Methods  summary 


Using  the  device  physics  models,  material  properties,  material 
costs,  manufacturing  costs,  and  heat  exchanger  costs,  the  metrics 
G  (Eq.  (1))  for  power  generation  and  H  (Eq.  (2))  for  thermoelectric 
cooling  can  be  evaluated  for  a  variety  of  traditional  and  novel 
thermoelectric  materials.  In  its  complete  form,  G  is  expressed 
as  [34] 
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where  the  junction  temperature  difference  is  evaluated  numeri¬ 
cally  by  solving  a  non-linear  system  of  coupled  equations  (Eqs. 
(3a)  and  (3b)).  A  point  of  diminishing  returns  in  system  cost  exists 
which  designates  the  optimum  fill  factor,  Fopt  [34].  The  true 
minimum  occurs  when  F  or  L  is  zero.  This  point  is  unrealistic  as 
no  device  would  exist.  Therefore,  the  diminishing  returns  point 
nearly  minimizes  G.  The  diminishing  returns  optimum  point  is 
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By  fixing  F  at  this  diminishing  returns  value,  Fopf,  an  optimum  L 
exists  that  uniquely  minimizes  G,  designated  as  Lop[.  We  present 
both  the  F=Fop[  and  F=1  scenarios  which  provide  the  optimized 
system  costs  at  an  optimum  LopC  given  a  fixed  F. 

The  analysis  provides  a  detailed  derivation  of  a  new  cost- 
performance  metric  G0 
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where  LT  is  the  ratio  of  the  material  thermal  conductivity,  k,  to 
the  heat  exchanger  thermal  conductance,  U.  G0  is  a  characteristic 
scaling  metric  for  G  and  accounts  for  material  cost,  application 
temperatures,  heat  exchanger  performance,  and  thermoelectric 
material  properties.  For  a  fixed  F,  when  the  metal  shunt,  ceramic 
plate,  heat  exchanger,  and  manufacturing  costs  are  excluded,  the 
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resulting  cost  metric  Gmin  is  the  minimum  material  cost  for  a 
thermoelectric  device  [34]: 

Gmin  =  4CqF  (14) 


The  fill  factor  can  equal  Fop[  or  some  other  value  as  determined  by 
device  design  requirements  not  considered  in  this  analysis.  This 
minimum  material  cost  scales  as  F2;  so  exceedingly  small  fill 
factors  will  allow  for  small  material  costs.  For  example,  for 
Material  #2  in  the  low  temperature  scenario  with  a  fill  factor  of 
F=0.01  the  minimum  material  cost  is  Gmin=$0.0084/W  which  is 
two  orders  of  magnitude  smaller  than  the  value  reported  in  Fig.  5 
with  Fopt=  0  .21.  This  minimum  represents  the  material  cost  domi¬ 
nated  regime;  however,  other  regimes  where  cost  is  dominated  by 
either  the  areal  manufacturing  or  the  heat  exchanger  also  exist 
[34]  and  have  different  minimums. 

For  cooling,  the  thermoelectric  operating  cost  H  is 


Ce  C"L2  +  C"L  J\+ZTm-\ 
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where  the  junction  temperature  difference  and  COP  are  deter¬ 
mined  by  solving  the  quadratic  in  Eq.  (7).  In  power  generation  and 
thermoelectric  cooling,  the  thermal  conductance  of  the  device,  the 
electrical  resistance,  and  the  overall  cost  are  complex  functions  of 
the  leg  length  L  and  fill  factor  F. 


3.  Results  and  discussion 

The  results  specifically  address  the  cost  of  candidate  materials, 
the  cost  to  process  those  materials  into  thermoelectric  legs  for  a 
module,  and  the  cost  of  the  heat  exchangers  making  up  the  full 
thermoelectric  device.  This  enables  realistic  assessments  of  cur¬ 
rent  and  future  materials'  potential  for  use  in  power  generation 
and  cooling  applications  by  linking  performance  with  estimated 
material,  manufacturing,  and  system  costs. 

3.1.  Material  and  module  manufacturing  costs 

The  contribution  of  manufacturing,  material,  and  heat  exchan¬ 
ger  costs  differ  by  the  manufacturing  and  material  types.  Fig.  4 
depicts  the  cost  for  processing  typical  bulk,  nanobulk,  nanowire, 


(Mat.  ID  #2)  (Mat.  ID  #3)  (Mat.  ID  #6)  (Mat.  ID  #14) 

Fig.  4.  Example  system  cost  breakdown  for  various  manufacturing  classes  on  a  $/W 
basis.  These  values  correspond  to  the  F — F opt  case  in  Fig.  5a.  The  low  temperature 
case  (Scenario  #1)  was  used,  and  the  cost  was  evaluated  for  the  optimum  fill  factor 
and  leg  length  (F=Fopt,  L=Lopt)  for  each  material.  For  nanowire  materials,  the  value 
$2500/m2  approximating  early  production  microfabrication  cost  was  used  in  the 
analysis  and  corresponds  to  about  $l/g  for  processing  silicon.  For  the  superlattice 
material,  $1800/m2  corresponding  to  $40/g  for  processing  lead  telluride  was  used. 
Areal  manufacturing  cost  (in  green)  represents  the  cost  of  dicing  and  the 
metallization  of  thermoelectric  material.  The  electrically  insulating  ceramic  plate 
and  heat  exchanger  costs  are  a  substantial  portion  of  the  overall  cost  for  all 
thermoelectric  devices.  (For  interpretation  of  the  references  to  color  in  this  figure 
caption,  the  reader  is  referred  to  the  web  version  of  this  paper.) 


and  superlattice  thermoelectric  materials  for  power  generation. 
The  numerical  values  and  cost  assumptions  are  presented  in  the 
Appendix  F.  A  similar  breakdown  can  be  performed  for  thermo¬ 
electric  cooling,  but  it  is  less  interesting  since  the  cost  of  electricity 
dominates.  For  power  generation,  the  cost  of  the  ceramic  plates 
and  heat  exchangers  are  a  dominant  cost  for  most  manufacturing 
classes.  Although  their  contribution  to  the  overall  cost  is  substan¬ 
tial,  heat  exchangers  can  increase  the  device  performance  since 
they  improve  heat  transfer  to  and  from  the  hot  and  cold  reservoirs. 
This  analysis  suggests  that  a  major  reduction  in  overall  system  cost 
can  be  achieved  by  making  inexpensive  heat  exchangers  (  <  $1.00/ 
(W/K))  with  high  U-values  (  >  200  W/m2  K).  Ceramic  plates,  typi¬ 
cally  used  as  substrates  in  planar  thermoelectric  modules,  provide 
the  necessary  electrical  insulation  and  mechanical  support  for  the 
device  and  represent  an  important  cost.  A  significant  reduction  in 
overall  system  cost  could  be  achieved  by  making  an  inexpensive 
(  <  $0.10/cm2),  thermally  conducting  (fc  >  40  W/m  I<),  electrically 
insulating,  and  structurally  suitable  substrate  material  that  can 
operate  at  the  given  application  temperature. 

The  advantage  of  nanobulk  materials  is  evident  in  Fig.  4  as 
illustrated  with  a  nanostructured  bismuth  telluride  material 
(Material  ID  #3)  compared  to  the  similar  bulk  material  (Material 
ID  #2).  The  nanobulk  material  has  a  lower  thermal  conductivity 
than  its  bulk  counterpart,  so  it  has  a  lower  $/W  value.  By  nanobulk 
processing,  a  hierarchy  of  length  scales  relevant  for  thermal 
and  electrical  transport  can  be  engineered  to  produce  a  better 
thermoelectric  material  [80].  These  processing  techniques  are  still 
affordable  and  demonstrate  the  benefit  of  a  slight  increase  in 
manufacturing  cost  that  results  in  an  overall  reduction  in  the  $/W 
value.  This  example  demonstrates  a  ~17%  overall  cost  reduction 
on  a  $/W  basis. 

Nanowire  materials  may  require  microfabrication  processes 
like  those  used  in  semiconductor  manufacturing.  In  semiconduc¬ 
tor  device  manufacturing,  cost  typically  scales  with  the  number  of 
fabrication  steps  involving  lithographic  patterning  and  etching,  so 
there  is  a  range  in  microfabrication  costs.  The  minimum  of  this 
range  corresponds  to  a  high-volume  process  similar  to  the  current 
state  of  silicon  photovoltaic  fabrication.  The  maximum  of  this 
range  corresponds  to  multi-mask  processes  used  to  make  com¬ 
plementary  metal  oxide  semiconductor  (CMOS)  devices.  In  the 
middle  cost  range  are  fabrication  steps  that  are  used  to  make 
micro-  electro-mechanical  systems  (MEMS);  this  is  the  cost  used 
for  nanowire  manufacturing  in  this  analysis. 

Superlattices  are  made  using  molecular  beam  epitaxy  (MBE) 
which  requires  many  hours  to  deposit  a  typical  thickness  of 
material,  so  the  cost  associated  with  achieving  sufficiently  thick 
thermoelectric  legs  is  high.  Since  the  manufacturing  costs  asso¬ 
ciated  with  thin  film  materials  are  significantly  higher  than  those 
for  bulk  materials,  the  increase  in  power  output  achieved  by  these 
materials  must  far  outweigh  the  additional  cost.  Given  the 
manufacturing  cost  of  MBE,  little  active  material  is  used  as 
demonstrated  by  the  small  leg  length  in  Table  1.  For  MBE  grown 
materials,  bulk  material  costs  are  low  because  little  active  material 
is  used.  The  cost  associated  with  material  deposited  on  surfaces 
other  than  the  active  substrate  is  not  considered  here  as  described 
earlier  in  Section  2.  Appendix  F  includes  details  regarding  a 
method  to  adjust  the  manufacturing  cost  and  compensate  for  a 
material  loss  factor. 

Metallization  and  dicing  costs  are  areal  manufacturing  costs; 
they  are  applied  to  all  materials  considered  herein  equally  since 
both  are  required  for  the  traditional  approach  of  assembling 
individual  legs  into  modules.  The  metallization  category  denotes 
the  deposition  of  diffusion  barriers  and  metallization  to  improve 
contact  to  the  electrical  shunt.  This  cost  includes  both  the  cost 
of  the  equipment  and  an  approximate  cost  of  the  metal  deposited. 
An  important  exception  to  traditional  metallization  and  dicing 
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processes  is  the  relatively  novel  polymer  thermoelectric  material 
processing.  The  metallization  and  formation  of  individual  leg 
structures  would  likely  be  done  entirely  through  screen  printing 
processes  or  a  similar  process  commensurate  with  the  unique 
fabrication  of  polymers.  The  impact  of  this  exception  is  evident  in 
the  final  G  values  computed  for  this  material  where  the  heat 
exchanger  cost  is  neglected,  and  the  system  parameters  used  are 
L=  100  mm,  F=  8,  TH=  100  °C.  The  screen  printing  approach  would 
lead  to  an  estimated  polymer  thermoelectric  module  cost  of  $140/ 
W  as  opposed  to  $550/W  if  traditional  metallization  and  dicing 
costs  were  applied. 

Large  mismatch  in  costs  can  yield  interesting  new  results. 
Specifically,  the  inexpensive  material  and  manufacturing  costs 
associated  with  the  polymer  and  its  low  thermal  conductivity, 
results  in  an  optimum  fill  factor  greater  than  one.  This  occurs 
because  the  cost  of  the  heat  exchanger  is  much  larger  than  the  cost 
of  the  material  (i.e.,  CHXU2t>C"'k').  While  a  fill  factor  greater  than 
one  does  not  adhere  to  traditional  thermoelectric  architecture,  it  is 
not  unimaginable.  This  is  simply  the  situation  where  the  cross 
sectional  area  of  the  legs  is  larger  than  the  plate  area.  This  is  the 
reverse  of  thermal  concentration  where  heat  is  actually  spread 
laterally  through  the  leg.  In  this  analysis,  the  same  heat  exchanger 
was  used  to  compare  all  materials.  This  heat  exchanger  (primarily 
its  cost)  is  not  suitable  for  polymers.  This  analysis  suggests 
polymers  require  different  geometries  and  heat  exchangers  than 
traditional  thermoelectrics. 

3.2.  Thermoelectric  power  generation 

The  cost-performance  analysis  for  thermoelectric  power  genera¬ 
tors  demonstrates  their  cost  competitiveness  when  using  different 
thermoelectric  materials.  Table  1  provides  an  overview  of  the 
material  types,  gravimetric  costs  ($/l<g),  properties,  and  optimal 
geometries.  In  this  table,  the  optimal  geometries  and  ZT  are  reported 
for  a  specific  temperature  application  (identified  in  Table  2  and 
defined  by  either  “low"  or  “high”)  where  that  material  performs  best. 
Other  temperatures  will  have  a  different  optimum  leg  length  and  fill 
factor  for  each  material.  In  some  cases,  particularly  nanowire  and 
superlattice  materials  ( e.g .,  Material  ID  #  6  and  14),  only  low 
temperature  material  property  data  is  available,  but  the  optimal  ZT 
values  may  be  at  higher  temperatures. 

There  are  complex  interplays  between  the  thermal  transport 
properties,  electrical  transport  properties,  and  costs  that  influence 
the  thermoelectric  leg  length.  For  example,  a  low  thermal  con¬ 
ductivity  material  allows  for  shorter  legs,  which  results  in  lower 
electrical  resistances  and  lower  volumetric  costs.  The  consequence 
of  high  thermal  conductivity  is  starkly  illustrated  by  nanobulk 
silicon  (Material  ID  #13)  where  the  high  thermal  conductivity 
(  >  10  W/m  K)  results  in  an  optimum  thermoelectric  leg  length  of 
over  90  mm.  A  high  electrical  conductivity  permits  longer  legs, 
which  results  in  larger  output  power  and  lower  $/W  cost.  There 
are  also  complex  interplays  between  manufacturing  areal  cost, 
heat  exchanger  cost,  heat  exchanger  thermal  conductance,  and 
material  thermal  conductivity  that  influence  the  fill  factor.  For 
example,  expensive  materials  favor  small  fill  factors  to  reduce 
costs  while  low  thermal  conductivity  materials  favor  larger  fill 
factors.  However,  the  fill  factor  is  commonly  fixed  by  other  system 
constraints,  and  an  optimal  thermoelectric  leg  length  is  sought. 

This  optimum  leg  length  is  the  length  that  minimizes  the  cost 
for  a  given  fill  factor  (typically  Fop[ )  and  is  determined  for  each 
material  and  temperature  scenario.  This  optimum  leg  length  is  not 
the  diminishing  returns  point  in  Eq.  (12)  but  rather  the  true  leg 
length  that  minimizes  the  system  cost  in  Eq.  (11)  for  a  fixed  fill 
factor.  The  optimum  leg  length  is  used  to  calculate  the  values 
reported  in  Figs.  4  and  5.  Three  fill  factor  cases  are  presented  in 
Fig.  5  for  each  of  the  temperature  scenarios.  The  first  case  is  for 


F=  1,  and  the  second  case  is  evaluated  at  the  optimum  fill  factor 
F=Fopt  (Eq.  (12)).  The  last  case  is  a  minimum  achievable  cost 
where  F  is  fixed  at  the  same  Fopt  value  as  the  second  case  (i.e., 
Gmm  =  4GoFop[2,  when  /<n  is  negligible).  This  minimum  cost  is 
reached  when  the  system  costs  of  ceramic  plates  and  heat 
exchangers  ( CHX )  as  well  as  the  metal  shunts  and  manufacturing 
costs  (C")  are  eliminated,  and  only  the  material  cost  (C'")  for  the 
thermoelectric  device  is  considered.  The  gap  between  this  mini¬ 
mum  cost  case  and  the  other  cases  represent  the  room  for 
improvement  in  manufacturing  and  system  costs. 

While  there  are  many  materials  that  can  operate  at  low 
temperatures  for  energy  scavenging  applications,  their  thermo¬ 
electric  materials  costs  alone  are  above  $1/W  (Fig.  5a).  The 
materials  with  the  lowest  minimum  cost  for  the  low  temperature 
scenario,  Tm  =  60  C,  are  predominately  bismuth  telluride-based 
chalcogenides  (Material  ID  #1-3)  although  the  bismuth  telluride 
nanowire  material  (Material  ID  #4)  has  high  costs.  The  presence  of 
certain  nanowire  (Material  ID  #14)  and  superlattice  (Material  ID 
#6)  materials  at  the  bottom  of  the  range  presented  for  the 
minimum  cost  in  Fig.  5a  highlights  the  potential  of  nanostructured 
materials  to  also  be  cost  competitive.  However,  the  gap  between 
the  minimum  cost  points  and  the  other  cases  indicates  that  the 
heat  exchanger,  ceramic  plate,  and  areal  manufacturing  costs  must 
be  reduced  (to  <  $1.00/(W/I<),  <  $1.00/cm2,  and  <$0.01 /cm2, 
respectively,  see  $ection  3.4)  for  most  thermoelectric  materials 
to  be  competitive.  One  option  to  reduce  manufacturing  costs  of 
nanomaterials  may  be  to  use  solution-synthesized  nanostructures 
since  they  do  not  require  specialized  microfabrication  equipment 
and  processes.  Nanostructured  bulk  material  made  with  processes 
like  ball  milling  also  have  lower  manufacturing  costs  than  nano¬ 
wire  and  superlattice  materials.  The  combined  effects  of  cost  and 
performance  are  responsible  for  this  result;  the  thermoelectric 
efficiency  improvements  outweigh  the  additional  costs. 

Certain  materials  are  clearly  promising  for  the  mid-  to  high- 
temperature  applications  as  shown  in  Fig.  5c.  The  lead  telluride- 
based,  half  Heusler,  skutterudite,  and  silicide  materials  designated 
by  Material  ID  #5,  7,  27,  21,  and  12,  respectively,  have  consistently 
low  costs  (  <  $6/W).  In  fact,  a  factor  of  two  reduction  in  the  system 
costs  for  these  materials  would  make  them  highly  competitive 
waste-heat  recovery  solutions  as  compared  to  Rankine  or  organic 
Rankine  cycles  at  $4-5/W,  especially  when  factors  such  as 
installation  costs,  reliability,  and  maintenance  are  also  considered. 
A  notably  high-cost  material  type  is  the  oxide  category  (Material 
ID  #23,  24).  In  spite  of  the  low  bulk  material  cost,  the  ZTm  is  too 
low,  largely  due  to  the  high  thermal  conductivity. 

These  results  should  be  placed  in  perspective  considering  the 
costs  of  competing  electricity  generation  technologies,  especially 
with  respect  to  the  heat  source  temperature  [7,81],  The  horizontal 
lines  in  Fig.  5  indicate  the  cost  of  competing  technologies.  When 
the  full  optimization  is  performed,  Fig.  5a-d  shows  that  current 
thermoelectrics  are  primarily  competitive  for  mid-  to  high- 
temperature  applications.  As  shown  in  Eq.  (13),  G0  is  related  to 
the  inverse  of  the  reservoir  temperature  difference  squared, 
( Th-Tc)~ 2,  which  makes  the  cost  of  thermoelectric  power  gen¬ 
eration  sensitive  to  the  reservoir  temperature  difference. 

Thermoelectric  power  generating  systems  are  more  costly  than 
primary  power  generation  sources  including  coal,  natural  gas, 
solar,  and  geothermal.  Primary  electricity  generating  sources 
such  as  coal  and  natural  gas  power  plants  cost  less  than  $3/W 
[82],  so  it  is  unlikely  that  thermoelectric  power  generation  would 
be  considered  as  an  alternative  to  these  power  sources,  especially 
considering  the  relatively  low  efficiencies  of  thermoelectrics 
compared  to  modern  Rankine  and  Brayton  cycles.  Through  cost 
reductions,  renewable  energy  technologies  have  recently  become 
more  competitive.  The  current  cost  target  for  photovoltaics  is 
'~$1/W,  and  the  photovoltaic  module  cost  is  projected  to  reach 
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Fig.  5.  Minimum  cost  of  thermoelectric  power  generation,  G,  in  $/W  for  materials  indicated  in  Table  1.  Comparisons  are  presented  for  the  four  scenarios  in  Table  2. 
Horizontal  lines  represent  costs  of  competitive  electricity  generation  technologies.  The  data  point  colors  represent  material  class  and  correspond  to  the  color  legend  in 
Figs.  6  and  7.  The  F=  1  line  represents  a  fill  factor  of  one;  the  F—Fopt  line  represents  a  module  design  with  the  optimum  fill  factors  from  Table  1.  At  these  same  F=F„pt  values, 
the  Material'  cost  line  (G=4G0FOpt2,  if  %  is  negligible)  represents  the  lowest  obtainable  cost  if  the  heat  exchanger  costs  and  areal  manufacturing  costs  are  considered 
negligible.  Equivalent  material-only  costs  for  other  fill  factors  can  be  estimated  using  the  Material's  cost  values  presented  here.  The  Fopt  values  tabulated  in  Table  1  can  be 
divided  out  to  yield  4G„  for  each  material,  and  another  value  for  F  can  be  inserted  14(T0Pj.  Hence  the  material-only  costs  can  be  smaller  than  those  plotted  by  more  than  an  order  of 
magnitude  for  the  smallest  practical  F  values  (~ 0.01  to  ~0.05).  (a)  Scenario  #1:  Th=  100  C,  Tc=  20  °C,  (b)  scenario  #2:  !«= 250  C,  Tc=20  °C,  (c)  scenario  #3:  TH= 500  C,  Tc=50  C  and 
(d)  scenario  #4:  TH=800  C,  Tc=50  °C.  (For  interpretation  of  the  references  to  color  in  this  figure  caption,  the  reader  is  referred  to  the  web  version  of  this  paper.) 


S0.50/W  [83,84],  Although  higher  at  S3-4/W,  the  costs  of  concen¬ 
trated  solar  power  [84]  and  geothermal  power  [82]  are  still 
competitive  and  preferable  to  the  costs  associated  with  thermo¬ 
electric  power  generation.  Key  exceptions  are  cases  where  com¬ 
pactness  and  portability  of  the  generator  are  critical,  and  many 
alternatives  to  thermoelectric  power  generation  are  not  feasible 
due  to  resource  availability,  size,  weight,  and  system  complexity. 

Thermoelectrics  are  especially  cost  competitive  compared  to 
other  waste-heat  recovery  technologies.  A  competitive  waste-heat 
recovery  technology  is  the  organic  Rankine  cycle  (ORC)  system  as 
shown  in  Fig.  5b.  The  cost  for  ORC  systems  is  approximately 
S4-5/W  [85],  Both  ORCs  and  thermoelectrics  can  be  used  to 
recover  unused  heat.  Thermoelectrics  can  even  be  complimentary 
to  ORCs  and  scavenge  additional  heat  not  recovered  by  the  organic 
Rankine  system  [86].  The  capital  costs,  system  footprint,  and 
maintenance  costs  for  organic  Rankine  systems  are  considerable, 
so  thermoelectric  devices  may  be  preferable  for  waste-heat 
recovery,  particularly  where  reliability  is  paramount. 

In  Fig.  6,  the  materials  can  be  compared  considering  two 
parameters  of  interest:  the  total  areal  costs  of  the  module: 

C"module=  C”'L0pt  +  C",  (16) 

and  the  figure  of  merit  ZTm.  The  dotted  reference  lines  are  drawn 
for  comparison  to  chalcogenide  materials:  bismuth  telluride  and 
lead  telluride  for  low  and  high  temperature  ranges,  respectively. 
These  plots  provide  an  image  of  the  current  materials  landscape 
and  can  illustrate  targets  for  future  materials  research.  There  is  a 
paucity  of  materials  in  the  lower  right-hand  quadrant  which 
indicates  new  material  system  modules  are  needed  that  are  less 
expensive  (C"modu/e<  $0.05/cm2)  and  higher  efficiency  (ZT >  1) 
than  the  chalcogenide  reference.  The  clustering  of  points  below 
the  chalcogenide  comparison  material  cost  point  indicates  sig¬ 
nificant  research  and  development  work  should  consider  material 


cost  while  improving  thermoelectric  efficiency.  For  example, 
nanowire-based  thermoelectric  materials  requiring  microfabrica¬ 
tion  processes  are  not  currently  competitive  options  as  shown  in 
Fig.  6b  and  c.  Materials  which  lie  in  the  upper  right-hand  quadrant 
may  be  desirable  in  some  instances;  military  applications  such  as 
mobile  power  sources  for  soldiers  or  space  applications  such  as 
space  vehicle  power  generators  are  two  examples. 

An  important  challenge  is  the  selection  of  materials  within  a 
given  quadrant.  When  faced  with  selecting  between  materials 
which  lie  in  relatively  close  proximity  to  one  another  in  Fig.  6, 
there  are  a  variety  of  application-specific  factors  that  govern  the 
selection  between,  for  example,  a  material  having  lower  cost  and 
performance  and  another  having  higher  cost  and  performance. 
Such  decision-making  will  likely  rely  on  additional  parameters 
related  to  feasibility  such  as  material  stability  with  respect  to  time, 
temperature,  and  compatibility  with  the  operating  environment, 
material  toxicity,  mechanical  properties,  and  compatibility  with 
metallization  and  other  module  materials,  none  of  which  are 
considered  in  this  analysis. 

3.3.  Thermoelectric  cooling 

Thermoelectric  cooling  is  an  application  where  thermoelectrics 
have  performed  well  in  the  market.  Much  of  this  success  can  be 
attributed  to  the  room  temperature  performance  of  bismuth 
telluride  alloys.  However,  other  thermoelectric  materials  may  also 
be  attractive  alternatives.  The  cost  analysis  discussed  previously, 
Eq.  (2)  and  (15),  was  performed  on  the  materials  listed  in  Table  1 
and  is  presented  in  Fig.  7.  Not  all  materials  are  able  to  achieve  the 
targeted  A7  =  l\  -  T2  associated  with  the  example  application. 
Eq.  (8)  shows  that  if  a  material  does  not  have  a  sufficiently  large 
ZT  for  a  given  AT,  then  the  COP  is  negative  indicating  that  cooling 
to  that  AT  cannot  be  achieved.  The  operating  cost  in  Fig.  7  is 
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expressed  for  each  material  at  the  optimum  leg  length  that 
minimizes  Eqs.  (2)  and  (15)  (reported  in  Table  F.2).  Materials  are 
listed  by  their  material  identification  number  and  are  ordered  by 
increasing  ZT.  There  is  a  clear  trend  that  larger  ZTs  have  lower 
operating  cost.  This  is  because  the  COP  is  the  strongest  driving 
parameter,  and  the  purchased  electricity  dominates  the  cost.  The 
coefficient  of  performance  is  directly  related  to  ZT. 

The  cost  of  thermoelectric  cooling  can  be  compared  against  the 
operating  cost  of  absorption  refrigeration  and  air  conditioning 
[14,87],  and  many  of  the  chalcogenides  are  cost  competitive.  Material 
ID  #6  has  the  potential  to  be  the  least  expensive  even  considering  its 
costly  MBE  process.  The  primary  challenge  is  making  high  quality 
superlattices  with  thickness  approaching  130  pm,  which  is  the 
optimum  length  as  determined  by  this  analysis.  Material  ID  #3  also 
shows  that  the  improvements  in  thermoelectric  power  output  and 
associated  cost  reduction  caused  by  nanostructuring  compensate  for 
the  additional  manufacturing  expenses.  Given  the  straightforward 
relationship  to  ZT,  an  ideal  material  with  ZTxo o  can  have  a  cost  of 
$0.004/l<W  hth  for  AT=10°C.  This  suggests  there  is  still  room  for 
material  improvements  in  cooling  applications  to  reduce  the  cost 
even  further.  Cost  reduction  can  be  achieved  by  increasing  ZT  with  a 
system  capital  cost  that  is  recovered  by  the  operating  cost  savings 
over  the  lifetime  of  the  device. 

3.4.  Thermoelectric  device  targets  summary 

The  key  targets  for  thermoelectric  devices  pertain  to  system 
parameters  and  costs.  Meeting  or  exceeding  these  targets  would  lead 
to  lower  cost  values  than  those  determined  here.  As  suggested  by  G0 
in  Eq.  (13),  both  large  ZT  and  small  volumetric  costs  C'"  are  desirable, 
but  they  need  not  occur  simultaneously.  Materials  with  low  ZT  and 
low  volumetric  costs  can  have  a  lower  $/W  value  than  those  having 
high  ZT  and  high  costs.  Hence,  there  is  no  unique  target  for  ZT  or 
volumetric  cost,  but  approximate  targets  can  be  set  for  current  state- 
of-the-art  system  components.  The  heat  exchanger  and  ceramic  plate 
dominate  the  overall  cost  for  many  devices.  For  the  situation  where 
the  heat  exchanger  costs  dominate,  improving  ZT  is  more  important 
than  reducing  C".  The  $1/(W/I<)  target  for  heat  exchanger  cost  is 
below  the  existing  minimum  cost  [68],  and  (T-values  above  200  W/ 
m2  I<  are  desirable.  Ceramic  plates  with  thermal  conductivity  above 
40  W/m  K  and  cost  below  $0.10/cm2  would  reduce  thermoelectric 
system  costs.  Since  the  areal  manufacturing  costs  contribute  mini¬ 
mally  to  the  overall  cost,  these  costs  should  remain  at  or  below  the 
current  value  of  $0.01 /cm2  for  dicing  and  metallization  of  bulk 
materials.  Any  additional  areal  processing  costs  would  increase  this 
value,  and  gains  in  other  parameters  would  need  to  mitigate  the 
additional  cost. 


4.  Conclusion 

This  work  applies  a  new  thermoelectrics  cost  analysis  [34] 
which  incorporates  material  properties,  device  physics,  material 
costs,  manufacturing  costs,  and  system  costs.  The  analysis  and 
resulting  cost  values  provide  a  tool  for  thermoelectric  device 
research  and  development.  Optimization  of  thermoelectric  leg 
length  and  fill  factor  minimizes  the  ratio  of  cost  to  performance, 
as  opposed  to  optimization  of  efficiency  and  power  output  alone. 
This  work  makes  progress  on  evaluating  the  status  of  thermo¬ 
electric  materials  relative  to  competitive  power  generation  tech¬ 
nologies  for  a  range  of  application  temperatures.  Novel  nanowire 
and  superlattice  materials  have  the  potential  to  have  a  low  $/W 
value  if  improvements  in  ZT  are  made  above  what  is  reported,  but 
with  the  currently  reported  values  they  are  not  competitive  in  the 
near-term  due  to  the  large  costs  associated  with  microfabrication/ 
MBE  manufacturing  techniques.  When  applied  to  several 


traditional  and  new  thermoelectric  material  classes,  the  analysis 
demonstrates  the  paucity  of  materials  which  are  both  higher 
performing  and  lower  cost  than  standard  chalcogenides.  However, 
there  are  multiple  thermoelectric  materials  which  could  yield 
power  generation  systems  that  are  affordable  and  competitive 
with  other  forms  of  power  generation;  the  key  challenges  for  these 
materials  will  be  in  both  engineering  devices  around  them  and 
scaling  their  production.  The  utility  of  the  cost  metrics  presented 
here  will  evolve  as  factors  such  as  resource  availability  and  energy 
costs  change  over  time  demonstrating  the  need  to  revisit  this 
analysis  periodically.  The  realization  of  commercial  thermoelectric 
devices  depends  on  the  simultaneous  improvement  and  optimiza¬ 
tion  of  material  properties,  system  design,  system  costs,  and 
material  and  manufacturing  costs. 
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